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Abstract 

First-order nonlinear interactions of Tol lmien-Schl ichting waves of 
different frequencies and initial amplitudes in boundary-layer flows are 
analyzed by using the method of multiple scales. For the case of two 
waves, a strong nonlinear interaction exists if one of the frequencies 
co 2 is twice the other frequency wi. Numerical results for flow past a 
flat plate show that this interaction mechanism is strongly destabilizing 
even in regions where either the fundamental or its harmonic is damped 
in the absence of the interaction. For the case of three waves, a 
strong nonlinear interaction exists when to 3 = 102- on. This combination 
resonance causes the amplitude of the wave with the difference frequency 
C03 to multiply many times in magnitude in a short distance even if it 
is damped in the absence of the interaction. The initial amplitudes 
play a dominant role in determining the changes in the amplitudes of 
the waves in both of these mechanisms. 
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I. INTRODUCTION 


One of the major roads from laminar to turbulent flow involves the 
initial linear amplification of disturbances, which might be present in 
the flow. However, as these disturbances grow to appreciable amplitudes, 
nonlinear effects set in. The nonlinear mechanisms that are activated 
depend on the spectrum of the disturbances. In this paper, we investigate 
two of these mechanisms. 

In his experiments on the transition from laminar to turbulent flow 

in a separated shear layer, Sato^ observed the appearance of the subharmonic 

of order one-half in addition to the higher harmonics of the fundamental 
2 

wave. Wille observed the development of subharmonic waves while investigating 

3 

the stability of both circular and plane jets. Kachanov, et al observed 

that, in addition to the higher harmonics of a fundamental wave, which 

was introduced in the flow by a vibrating ribbon, a subharmonic wave 

with one-half the frequency of the fundamental wave appeared downstream. 

4 

Michalke postulated that the subharmonic appears when two vortices 

5 

rotate around each other in a fusion mating dance. Kelly showed that 
the appearance of the subharmonic in a shear layer is due to a secondary 
linear instability associated with a time-dependent flow that consists 
of the superposition of the basic flow and a finite-amplitude funda- 
mental wave. Nayfeh and Bozatli^ investigated the appearance of the 
subharmonic in boundary layers by analyzing the instability associated 
with a time-dependent flow that consists of the superposition of the 
basic flow and a Tol Imien-Schl ichting wave. The results show that the 
amplitude of the fundamental wave must exceed a critical value to 
trigger this parametric instability. This value is proportional to a 
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detuning parameter that is the real part of k - 2K, where k and K are 

the wavenumbers of the fundamental and its subharmonic, respectively. 

For the Blasius flow, the critical amplitude is approximately 29% of the 

mean flow. For other flows where the detuning parameter is small, such 

as free-shear layer flows, the critical amplitude can be small, thus the 

parametric instability might play a greater role. Since the analysis of 
5 6 

Kelly and Nayfeh and Bozatli are linear, they do not account for the 

effect of the subharmonic wave on the fundamental wave. This effect may 

be small initially, but as the subharmonic grows appreciably, its 

effect on the fundamental cannot be neglected. One of the purposes of 

the present paper is to determine the nonlinear interaction of a Toll mi en- 

Schl ichting wave with its subharmonic. 

7 8 9 

Sato , Miksad , and Kachanov, et al observed that the nonlinear 
development of the waves in the transition region depends on the initial 
and external disturbances. Sato 7 conducted an experiment on the stability 
of symmetric laminar wakes by exciting two unstable modes with the 
frequencies fi and f 2 . He observed the generation of waves having the 

O 

frequencies f 2 ± fi. Miksad excited two unstable modes of a laminar 
asymmetric free-shear layer. He also observed nonlinear triggered 
instabilities of the difference mode f 2 - fi, subharmonics, and higher 
harmonics of the fundamental waves. Kachanov, et al introduced two 
Tollmien-Schl ichting waves in the boundary layer on a flat plate by 
using two vibrating ribbons. They observed the appearance and growth of 
a Tol Imien-Schl ichting wave having the difference frequency f 2 - fi. 
Norman^ also observed the amplification of the difference harmonic of 
two introduced disturbance waves in his experimental study of secondary 
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flows around and downstream of protuberances in laminar boundary layers. 
The second purpose of the present paper is to determine the nonlinear 
interaction of three Tol lmien-Schl ichting waves (combination resonance) 
in boundary layers and show that the difference frequency can be very 
unstable when generated by the nonlinearity, even though it is stable 
when introduced by itself in the boundary layer. 

The problem is formulated in Sec. II. The analysis for the com- 
bination-resonance case is contained in Sec. Ill, while the results for 
the second-harmonic case are stated in Sec. IV. The numerical procedure 
is discussed in Sec. V, while the numerical results are presented in 
Sec. VI. 
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II. PROBLEM FORMULATION 


We consider nonlinear interactions of wave packets in a two-dimensional 
steady incompressible boundary-layer. The equations describing the 
motion of the fluid are 


3u 
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where 


,2 = If. 

ax' 


ay' 


Here, x and y are made dimensionless by using a reference length 6 r , 
the time is made dimensionless by using dy/U^, and the velocities are 
made dimensionless by using the freestream velocity U^. The Reynolds 
number R = UyJy/v with v being the fluid kinematic viscosity. 

The analysis is restricted to basic flows that are slightly 
nonparallel (i.e., vary slowly in the streamwise direction) and to 
disturbances that are small but finite. The slow variation is expressed 
by using the slow scale Xi =eix, where ei is a small dimensionless quantity 
that characterizes the nonparallelism of the flow and can be related 
to R by ei = R" 1 . The smallness of the amplitude of the disturbance 
is expressed by introducing the small dimensionless parameter e. For 
a general solution, we assume that e = 0(ci) so that the resulting 
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expansion accounts simultaneously for the effects of nonparallelism and 
nonlinearity. When e « ei, the nonlinear effects are negligible and 
the solution reduces to those obtained in Refs. 11 and 12. When 
e » ei, the nonparallel effects are negligible and the solution reduces 
to equations with constant coefficients. 

We assume that each flow quantity is the sum of a mean-flow quan- 
tity and an unsteady disturbance quantity, which is assumed to be much 
smaller than the mean-flow quantity. We can then express the velocity 
components and the pressure as 

u(x,y,t) = U 0 (xi,y) + eu(x,y,t), (6) 

v(x,y,t) = eiV 0 (xi,y) + ev(x,y,t), (7) 

p(x,y,t) = P ( Xl ) + ep(x,y,t), (8) 


where U 0 , V 0 , and P 0 are the nonparallel basic-flow quantities. Sub- 
stituting Eqs. (6) - (8) into Eqs. (1) - (5) and subtracting the basic- 
flow quantities, we obtain 
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Without loss of generality, we let e = ei. To determine the wave- 

packet solutions of Eqs. (9) - (13), we use the method of multiple 
13 

scales and seek an expansion in the form 


u = Uo(x 0 ,Xi,y,T 0 ,Ti) + 

eu^Xo.Xi.y.ToJi) + .... 

(14) 

v = Vo(x 0 ,Xi,y,T 0 ,Ti) + 

eVi(x 0 ,Xi,y,T 0 ,Ti) + .. ., 

(15) 

p = Po(x 0 ,Xi,y,T 0 ,Ti) + 

£-Pi ( xo » x i ,y , T o , T i ) + ..., 

(16) 


where xo = x, T 0 = t, and Ti = et. Substituting Eqs. (14) - (16) into 
Eqs. (9) - (13) and equating coefficients of like powers of e, we obtain 
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where 



In what follows, we describe the details of the analysis for the 
combination-resonance case and only state the results for the second- 
haromnic resonance case. 
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III. COMBINATION RESONANCES 


A. First-Order Problem 

For the case of combination resonances, we consider three wave- 
packets centered at the frequencies 101 , 012 , and co 3 , Then we examine the 
resonances that might exist among them. Thus, the solution of Eqs. 

(17) - (21) is expressed as a linear combination of three Tollmien- 
Schlichting waves; that is,. 

u 0 = A 1 (x 1 ,T 1 )Cn(y;xi)exp(i0i) + A 2 (xi,Ti) x 

Ci2(y;xi)exp(i0 2 ) + A 3 (xi ,T i)ci 3 (y;xi)exp(i0 3 ) + c.c., (27) 
v 0 = Ai (xi ,Ti )C 2 i (y*»Xi )exp(i0i) + A 2 (xi,Ti) x 

C22(y;xi)exp(i0 2 ) + A 3 (xi,Ti)c 23 (y;xi)exp(ie 3 ) + c.c., (28) 

Po= A 1 (xi,Ti)c 31 (y;xi)exp(i9i) + A 2 (xi,Ti) x 

C32(y;xi)exp(i9 2 ) + A 3 (xi,Ti)t 33 (y;xi)exp(i0 3 ) +c.c., (29) 

where 

30 30 

3x7 = k n^ Xl ^ ’ 7T7 = " % ( n = 2 » 3 ) (3°) 

with the oo n being real constants. The quasi -parallel Orr-Sommerfeld 
problems for these waves are 

M 'Kn- t 2 „i k n> = D \o * iklt ,n = °- (3,) 

C 2 „- V k n' “n> = i < u » k „ * + ? 2 n 0U » 

* ik -? 3 n - 1 <° 2 - k n>S„ * »• ‘ 32 > 
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where D = 8/3y. 

B. Second-Order Problem 

Substituting Eqs. (27) - (29) into Eqs. (22) - (26), we find that 
the inhomogeneous parts in Eqs. (22) - (26) contain terms proportional 
to 

exp(i 0 i), exp(ie 2 ) , exp(i© 3 ), 

exp[i (0 2 - Fi)], exp[i(0i + 0 3 )], exp[i(0 2 - 'e 3 )] 

where the overbar indicates the complex conjugate. The terms that are 
proportional to these exponential expressions will create secular terms 
in the particular solutions for ui, Vi, and p 3 if k 3 T k 2 - k 3 
and to 3 ~ o) 2 - tox ; that is, when a combination resonance exists among the 
waves. To express quantitively the nearness of the above resonances, we 
introduce the two detuning parameters and a 2 defined by 

oj 3 - w 2 +(*)!= ecri , (36) 

Real(k 3 - k 2 + k 3 ) = ea 2 , (37) 
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where a 

n 

= 0(1 ). Using Eqs 

. (36) 

and (37), 

we write 
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where k^ stands for the imaginary part of and 

<f> = J" a 2 dx x - aiTi (41 ) 

To determine the A n , we seek a particular solution for the second- 
order problem in the form 

ui = ihi(y;xi)exp(i6i) + ^i 2 (y;xi)eyp(ie 2 ) +i|>i 3 (y;xi) x 

exp(i0 3 ) + c.c. , (42) 

Vi = ^ 2 i(y;xi)exp(i 0 i) + ^ 22 (y;xi)exp(i0 2 ) + i|> 23 (y;xi) x 


exp(i0 3 ) + c.c. , (43) 

Pi = 'bi(y;xi)exp(i0i) + ip 3 2 (y ; xi )exp(i0 2 ) + ip 3 3 (y ;xi ) x 

exp(i0 3 ) + c.c. (44) 


Substituting Eqs. (27) - (30) and (38) - (44) into Eqs. (22) - (26) and 
equating the coefficents of exp(i0j), exp(i© 2 ), and exp(i0 3 ) on both 
sides, we obtain the following equations: 
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for j = 1, 2, and 3, where the d. . are given in Appendix A. 

* J 
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C. Adjoint Problem 

Since the homogeneous parts of Eqs. (45) - (49) are the same as 
Eqs. (31) - (35) and since the latter have a nontrivial solution, the 
inhomogeneous equations (45) - (49) have a solution if, and only if, the 
inhomogeneous' parts are orthogonal to every solution of the adjoint 
homogeneous problem; that is. 
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we obtain the following differential equations for the evolution 
of ai , a 2 , and a 3 : 


J- 1 . 8a i _ 
to i 8t 3x 


J_ 8 a 2 , 8 a 2 
gj 2 8 t 8 x 


1 _ 8a 3 + 833 _ 

(O 3 8 t 8 x 


(ei - k^Jai + 

ei -Sr 2 - - k . )a 2 + 
T 2 2 * 

£1 •pj 2 ■ k^. )a 3 + 


e a 2 a 3 exp(-i<j>), 

e ir 211 - aia 3 exp(1<t)), 

1 2 

£ S 222 a 2 a’iexp(-i(j)) , 
1 3 


(57) 

(58) 

(59) 


where w' = du) n /dk n is the group velocity and <j> is defined in Eq. (41). For 
spatial modulation only, ai = 0 and 8 a n / 8 t = 0 ; all the calculations 
presented in this paper are for this case. We note that Eqs. (57) - 
(59) account for the combined effects of the nonparallelism (i.e., 
growth of the boundary layer) and the nonlinear interaction. If 
e « ei, the nonlinear interactions can be neglected and the spatial 
variations in Eqs. (57) - (59) reduce to the nonparallel solutions of 
Refs. 11 and 12. When ei « e, the effects of the nonparallelism are 
negligible; that is, one can set ei = 0 and all the coefficients in Eqs. 
(57) - (59) can be treated as constants. 


13 



IV. HARMONIC RESONANCE 


The interaction between a fundamental Tollmien-Schlichting wave and 
its second harmonic is analyzed using a procedure similar to that 
outlined in the previous section. In this case, instead of Eqs. (57) - 
(59) we obtain 


1 4. 1 

oji 9t 9x 


k .)ai 

lT 


(60) 


1 932 , 932 

0)2 9t 9x 





a?exp(i<J>) , 


(61) 


where <f> is defined in Eq. (41) and 


eo 2 = Real(k 2 - 2ki), eai= 0)2 - 2u>i , (62) 

and fi, f 2 , hn, h 2 2 > hi 2 , and h 2 i are given in Appendices B and C. 

For spatial modulation only, Oi = 0 and 9a n /9t = 0. All the 
calculations presented in this paper are for this case. 
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V. COMPUTATION PROCEDURE 


A. Solutions of First- and Second-Order Problems 

The same procedure is followed in solving the first- and second- 
order problems for both harmonic and combination resonances. Therefore, 
only the computation methodology for the solution of the first-order 
problem for the first mode is outlined here. 

Equations (31) - (33) are expressed as a system of first-order 
differential equation in the form 

at - Gz - < 63 > 

where z is a 4 x 1 matrix with the elements 

zi = Cu(y;xi), z 2 = Dei i (y;x x ) , z 3 =4 2 i(y;xi), z 4 = c 3 i(y;xi), 

(64) 

and G is a 4 x 4 matrix; its elements are given in Appendix D. 

We start the integration of Eqs. (63) at y = y , where y £ is larger 
than the boundary-layer thickness. Hence, U 0 = 1, DU 0 = 0, and 
D 2 U 0 = 0 at y . Then the matrix G has constant coefficients at y = y g 
and Eqs. (63) have solutions of the form 

4 

z, = l c--exp(A.y) for i = 1, 2, 3, and 4, (65) 

3 -i — * J J 

J-l 

where the c.. are constants, the A's are the solutions of 

* *J 

|G - AI| = 0, (66) 

and I is the identity matrix. Equation (66) has the roots 

A = ± k i , A = ± [k 2 + i ( k x - mi)R] ^ 2 , (67) 

1 >2 3 > 4 1 
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Two of these roots have positive real parts that make the solution grow 
exponentially as y -*• hence, they must be discarded to satisfy conditions 
(35). This leaves two linearly independent solutions that decay ex- 
ponentially with y. 

The eigenvalues are not known a priori and must be determined along 

with the eigenfunctions. For given values of wi and R, we guess a value 

for ki and integrate Eqs. (63) from y g to y = 0. If the quessed value 

of ki does not staisfy the boundary conditions at y = 0, k! is incremented 

by using a Newton-Raphson scheme and the procedure is repeated until the 

boundary conditions are satisfied to within a specified accuracy. The 

integration is done by using a computer code developed by Scott and 
14 

Watts . This technique orthonormal izes the solution of the set of 
equations whenever a loss of independence is detected. 

B. Solution of Adjoint Problem 

The solution procedure is exactly the same as that for the first- 
order problem. The coefficients of the z matrix are 

ic k k 

zi = C 2 i(y;xi), z 2 = Dc 2 i(y;xi), z 3 = 4 3 i(y;xi), 

k 

z 4 = Cn(y;xi) (68) 

and the adjoint problem has the same eigenvalues as the first-order 
problem. 

C. Solvability Conditions 

The calculations are repeated at different streamwise locations to 

evaluate f., h-., k., and the other interaction integrals for a given 
J J J J 

frequency along the x-axis. A fourth-order fixed step-size Runge-Kutta 
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integration scheme is used to solve either Ecs. (57) - 
bination resonances or Eqs. (60) and (61) for harmonic 
find the amplitudes of the waves for different initial 
respective modes. 


(59) for corn- 
resonances to 
amplitudes of the 


17 



VI. Results and Discussion 


The analysis presented in this paper is applicable to both two-and 
three-wave interactions. First, we present and discuss numerical 
results for the case of two-wave interactions. Then, we present and 
discuss numerical results for the interaction of three waves whose 
frequencies are such that F 3 = F 2 - F x . 

A. Two-Wave Interactions 

The numerical results presented in Ref. 6 show that the amplitude 
of a wave af = ea 2 must exceed a critical value before it can generate 
and amplify its subharmonic. For the Blasius flow, the critical value 
is approximately 29% of the mean flow. This is for the case when the 
subharmonic wave has an infinitesimal amplitude. When the amplitude 

k 

ai = eai of the subharmonic wave is not infinitesimal, its influence on 

k 

a 2 should be taken into account. The equations governing this influence 

are Eqs. (60) and (61) whose general solution is not available yet. The 

6 * 

previous results of the parametric instability model show that a x 

k 

oscillates about its non-interaction value until a 2 reaches the critical 
value. Figures 1 and 2, obtained by numerically solving Eqs. (60) and 

k k 

(61), agree with this conclusion. Initially, a 2 increases while ai 

oscillates around its non-interaction value. 

★ 

At R ~ 580, Fig. 1 shows that ai starts to deviate sharply from its 

* 

non-interaction value, while it follows from Fig. 2 that £na 2 ~ -1.25 

k k 

or a 2 = 0.286 at this location. Hence, when a 2 is less than this 
★ 

critical value, ai can be approximated by its non-interaction value; 
that is 
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( 69 ) 


* * 

ai = ai 0 exp(-k il .x + it) 

★ 

where ai 0 and x are the initial amplitude and phase of the subharmonic 

wave, respectively. If we substitute Eq. (69) into Eq. (61) and neglect 

the nonparallel effects, we obtain 
★ 

+ k zi a* = at§exp[-(2k i1 + ica 2 )x + 2it] - (70) 

The solution of Eq. (70) that satisfies the initial condition 
* * 

a 2 = a 2 (0) at x = 0 can be written as 

a* = [a*(0) + f j 2k h | 1 ieg2 ) at§exp(2i T )]exp(-k 2 .x) 

1 i 

- f ; T 2k fo " eoV) ~ a ^exp[-(2k ii + ieo 2 )x + 2ix]. 

1 1 

(7!) 

Equation (71) represents an approximation to a 2 as long as it is less 
than the critical value needed to trigger the parametric instability in 
the subharmonic wave. 

Next, we consider the generation and amplification of a second- 
harmonic wave by a fundamental Tollmien-Schlichting wave. We consider 
the following three cases: (i) fundamental wave is stable while its 

second harmonic is unstable, (ii) fundamental wave is unstable while its 
second harmonic is stable, (iii) both fundamental and second-harmonic 
waves are unstable. 

When the fundamental wave is initially stable while its second 
★ 

harmonic is unstable, ai decays until it reaches the unstable region and 

then it increases as shown in Fig. 1. For Reynolds Number less than 
★ 

560, ai oscillates around its non-interaction value, implying a small 

★ 

initial influence of its second harmonic on it. Thus, a 2 can be approx- 
imated initially by Eq. (71). Figure 2 shows that the values obtained 
from Eq. (71) are in good agreement with those obtained by numerically 
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integrating Eqs. (60) and (61) for R _< 560. After a short initial 

distance, the second term on the right-hand side of Eq. (71) decays. 
★ 

Then, a 2 can be approximated by 


★ 

a2= 



( 0 ) + 


h 2 1 

f 2 (2k ii + i ea 2 ) 


a 10 exp(2iT)]exp(-k .x), 

2 » 

(72) 


as long as a 2 is less than the critical value. 

Hence, the effect of the fundamental wave on its second harmonic is 

ic 

to increase its initial amplitude. However, as a 2 attains large values, 

ic ic 

it strongly influences ai which in turn strongly influences a 2 . The 
result is an accelerated instability. 

For the case when the fundamental wave is initially unstable while 
its second harmonic is stable, we performed calculations for waves with 
the frequencies Fi = 46.5 x 10" 6 and F 2 = 93 X 10 -6 . The fundamental 
wave is in the unstable region at R = 950 where the calculations are 
started. Thus, its unstable downstream of R = 950. Figure 3, obtained 

ic 

by numerically integrating Eqs. (60) and (61), shows that ai hardly 

★ 

deviates from its non-interaction value. On the other hand, a 2 in- 
creases many orders of magnitude even for small initial amplitudes of 
the fundamental wave as shown in Fig. 4. In these calculations, the 
initial amplitude of the second-harmonic wave is taken to be 0.1% while 
the initial amplitudes of the fundamental wave are 0.1% and 0.5%. Since 

ic 

ai hardly deviates from its non-interaction value, Eq. (71) is expected 

★ 

to be a good approximation to a 2 . Figure 4 shows that the values obtained 
from Eq. (71) oscillate about those obtained by numerically integrating 

ic 

Eqs. (60) and (61). Since the initial values are very small, a 2 

★ 

does not reach the critical value to influence ai. After a short initial 
distance, the first term on the right-hand side of Eq. (71) decays and 
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★ 

a 2 can be approximated by 

"k 2 

32 = " f 2 (2k - /-Mega) aioexp[-(2k • + iea 2 )x + 2ix], 

i i 

( 73 ) 

★ 

as long as a 2 is less than the critical value. Consequently, the effect 
of the interaction is to produce a second-harmonic wave that grows 
approximately at a rate that is twice that of the fundamental wave. 

For the case when both waves are unstable, we performed numerical 
calculations for waves having the frequencies Fi = 52 x 10”^ and 
F 2 = 104 x 10"^ starting near R ~ 600. Figure 5 shows that initially 
aj deviates slightly from its non-interaction values. Hence, Eq. (71) is 

■k 

expected to be initially a good approximation to a 2 . Figure 6 shows 

that the numerical values obtained from Eq. (71) are in good agreement 

with those obtained by numerically integrating Eqs. (60) and (61) when 
★ 

a 2 is less than 0.29. Thus, in this case, the effect of the interaction 
on the second-harmonic wave is to increase its initial amplitude and to 
produce a term that grows at a rate that is twice the growth rate of 
the fundamental wave. Due to the fact that both waves are initially 
unstable, the interaction is more effective in this case than in the 
preceding two cases. 

B. Three-Wave Interactions 

In their experimental studies, Kachanov, et al^, Miksad®, Norman^, 
and Sato' 7 introduced two separate waves of different frequencies into 
the flow that was being studied. They observed the growth of a wave 
whose frequency is equal to the difference frequency. Kachanov et al 

-fi -fi 

used the frequency pairs Fi = 88 x 10” D and F 2 = 104 x 10 and 
-fi -fi 

Fi = 88 x 10 and F 2 = 120 x 10 to analyze the growth of the associated 
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difference-harmonic waves at F 3 = 16 x 10 “ 6 (i.e. F 3 = F 2 - Fi) and 
F 3 = 32 x 10"^ in a boundary-layer flow over a flat plate. Using the 
same frequency pairs, we determined the amplitudes of the fundamental 

ic "k k 

waves (a 3 and a 2 ) and the difference-harmonic wave a 3 by numerically 
solving Eqs. (57) - (59). Figures 7 and 8 show the large increases in 
the amplitudes of the above difference-harmonic waves due to the nonlinear 
interaction of the waves with the frequencies Fi and F 2 . The amplitude 
of the difference-harmonic wave at R ~ 715 is increased 120-200 times 
its non-interaction value. 

However, the amplitudes of the fundamental waves change very little 

from their non-interaction values as shown in Fig. 9 for F 3 = 88 x 10”^. 

★ ★ 

Thus, ai and a 2 can be approximated by 

•k -k kk 

ai = ai 0 exp(-k .x + in), a 2 = a 20 exp(-k -x + ix 2 ), (74) 

11 2 1 

* * 

where am and a 20 are the initial amplitudes and xi and t 2 are the 

initial phases of the fundamental waves. Substituting Eq. (74) into 

Eq. (59) and neglecting the nonparallel effects, we obtain 
★ 

^ + k^ .a* = at 0 a 20 exp[-(k ji + k^. + iea 2 )x 

+ i (ti + x 2 )] (75) 

The solution of Eq. (75) that satisfies the initial condition 

k k 

a 3 = a 3 ( 0 ) at x = 0 can be expressed as 

*k 'k U k k 

a 3 = [a 3 (0) + f ( k . ^ ; + -” e g2 ) ai 0 a 20 exp(iTi + ix 2 )] x 

i i 2 i 

exp(-k Sl x) - f (k h »\ - 4 . , ) at.at,exp[-(k , + k ( 

1 ' 2 I 

+ iea 2 )x + i(n + x 2 )]• (76) 

k 

Equation (76) represents an approximation to a 3 as long as the amplitudes 
of the fundamental waves do not deviate from their non-interaction values. 
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Figures 7 and 8 show that Eq. (76) is initially in good agreement with 
those obtained by numerically integrating Eqs. (57) - (59). 

According to Eq. (76), the difference-harmonic wave grows at a rate 
that is the sum of the growth rates of the fundamental waves. Since the 
fundamental waves are unstable at R = 430, where the calculations are 
started, the difference-harmonic wave amplifies considerably, inspite of 
the fact that it is stable in the absence of the interaction. These 
results are in qualitative agreement with the experimental observations 
of Refs. 7 - 10. 
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APPENDIX A 


9A. 3?i .• 

dl j = ( 9x7 Cl J + A j 87T )j 


d2j “ ^2jO + d 2 j 1 » 


d2jO“ “(UoClj + ?3j 


3?3 A 


1 + |i(k 


O,* 3A. 9d, 

R“ k j Cl j ‘ VoD?1 j ) 97?- + [ ‘ Uo W 1 


3 ?i ,* dk. 

— 1 


9xx R v j 9Xi dxi j 9Xi ^j J j : 

d 2 ii = [i (k 2 - k 3 )ci 2 ?i 3 + C22DC13 + DC12C233A2A3 x 

exp[-i<f> ~ j + k 3 -j - ki-)dx], 

d 221 = C i ( k 1 + k 3)?ll?13 + ? 2 lD^i 3 + Dc 1 1 C 2 3 ] A 1 A 3 X 

exp[i<j) - / (k Xi + k 3i - k 2i )dx], 
d231 = ( - i k i £ I 1 C 1 2 + ik 2 Cl2?l 1 + ?22D?11 + D^ 12 ^2l)^2^1 X 


exp[- i 4> -J[ kx- + k 2i - k 3i )dx]. 


d3 j = d 3 j° + S 1 ’ 

2i 


9A. 


9^2 ■ 


d3 j 0= - (U °^j - IT k j^j } 3^ - [U ° 9^ + VoD?2 J + ?2 J DVo 

dk, 9? 2 , 

" R~ ^ dxT ^ 2 j + k j 9Xj ^ A j’ 


1 3 1 1 


*321 


(-ik 3 ?l 2 ?23 + ik 2 ? 22?13 + 522^^23 + ^ 22 ^ 23 )^ 2^3 X 

ex P[-i 0 -/(k 2i + k 3 . - kj i )dx] , 

(ik 3 ?i x? 23 + ik 1 ^ 2 iCi 3 + ^ 21^23 + D^ 2 iC 23 )AiA 3 x 

exp[i<f) - yjki- + k 3i - k 2i )dx]. 


(Al) 

(A2) 


(A3) 


(A4) 


(A5) 


(A6) 

(A7) 


(A8) 


(A9) 


(A10) 
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{ -ikiCzi + "i k 2 C2 2 C 1 1 + C22DC21 


+ 0^22^21)^2^1 



APPENDIX B 


f j - 1 (- - 


2ik . * 2ik . 

[ (Uo p )Cn + C3,-]C2< - (Uo - -D- 1 ) 


R 


C 2j C 3j dy 


( j = 1, 2, 3), 


JJ 


oo / 

■IP 


dk. 


Zii + [U 0 


+ «>j) + fr + - ir 1 > ht 1 +V » D ^ 


9£i,- 9^3. 


2i 


3Xi 9Xi 


+ Vo D?n - ~ (k. 


9?i . 


R 3Xi 


2ik . 3?2 ,• 


(Bl) 


2i ^-i 

+ ^ DV ° ■ % ^ 


' 

Ca Jct-j 

1 J J 


dy (j = 1, 2, 3), 


1123 


71 


__ ^ ^ 

[i(k 2 - k3)Cl 2 Cl3 + £ 22 DCl3 + DCl 2 C 2 3]C 2 l 


+ ( “ik3£l2C23 + ik2^ 22 ?13 + ^ 22 D^23 + D^2 2 4 2 3 ) ^ 3 1 


1 2 1 3 


7 ! 


[i(ki + k 3 )Cll?13 + ^2lDCl 3 + DCn?2 3 ]? 22 


+ {ik 3 ?l 1^23 + iklC 2 lCl3 + C 2 lD?23 + DC 2 lC 2 3)?32 


>312 


n 


| dy, 

| dy, 


[i (k 2 - ki )Cl2?ll + ? 22 D?11 + D^i 2 C2l]?23 


(B2) 


(B3) 


(B4) 


+ (-iki£ 2 iCi2 + ik 2 C22Cn + C 22 D^2i + D^22^2i)?3 3 !dy . (B5) 
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APPENDIX C 



+ [ 1 l<2 C2 2 ^ 1 1 

OO 

[(IkiCi! 



__ -jlf 

+ C21DC12 + i(k 2 - 

ikVS21?12+ ?22Dl 2 l + ^2lD?2 2 ]C31 | 

+ C2lDCn)?J 2 + (lkiCnC2i + C2lDC2l)C 



APPENDIX D 


g 1 1 = 0 , gi2 = 1, gi3 = 0, gi4 - 0 (Dl) 
g 2 i = i(U 0 k - u)R + k 2 , g 22 = 0, g 23 = R y 1 , g 2 4 = ikR (D2) 
931 = - i k, g 32 = g 33 = g 3 4 = 0 (D3) 
941 = 0, g 42 = -ik/R, g 43 = - [ i ( U 0 k - u) + k 2 /R], g 4 4 = 0 (D4) 

Q is a 2 x 4 matrix consisting of the last two rows of the matrix 


The matrix B has the 

elements: 




bn 

= bi 2 = b 13 = bn 

= 1 




(D5) 

b 2 1 

= -k, b 22 = k. 

b 23 = 

k, b 2 4 = - 

k 


(D6) 

b 3 1 

= i, b 32 = ik/k. 

b 33 

= -i, b 34 

= . 

-ik/k 

(D7) 

b4 1 

= (to/ k - 1), b 42 

= 0 , 

b4 3 = (<Vk - 

1): 

) b44 = 0 

(D8) 


where 

k = [k 2 + i(k - to) R3 1/2 (D9) 
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Amplitude of a wave at Fi = 52 x 10” 6 involved 
in a subharmonic resonance with a fundamental wave 
at F 2 = 104 x 10" 6 . 


Figure 1. 



In (a?) 


t 



Amplitude of the fundamental wave at F 2 = 104 x 10"® 
involved in resonance with its subharmonic at F x = 52 x 10“ 6 . 


Figure 2. 
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Figure 3. Amplitude of a fundamental wave at Fi = 46.5 x 10~® 

involved in resonance with its harmonic at F 2 = 93 x 10" 6 . 
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Figure 5. Amplitude of fundamental wave at Fi = 52 x 10 ® 
in resonance with its harmonic at F 2 = 104 x 10 


775 


involved 
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Amplitude of the second-harmonic wave at F 2 = 104 x 10 - ^ 
when involved in resonance with a wave at Fi = 52 x 10"^. 


Figure 6. 
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generated from a combination resonance of Fi = 88 x 10" 6 
and F 2 = 104 x 10" 6 . 
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involved in a combination resonance with two waves 
at F, = 104 x 10' 6 and F a = 16 x 10" 6 . 
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